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ABSTRACT 

We study the applicability of several galaxy environment measures (n th -nearest-neighbor distance, 
counts in an aperture, and Voronoi volume) within deep redshift surveys. Mock galaxy catalogs are 
employed to mimic representative photometric and spectroscopic surveys at high redshift {z ~ 1). We 
investigate the effects of survey edges, redshift precision, redshift-space distortions, and target selection 
upon each environment measure. We find that even optimistic photometric redshift errors (a z = 0.02) 
smear out the line-of-sight galaxy distribution irretrievably on small scales; this significantly limits the 
application of photometric redshift surveys to environment studies. Edges and holes in a survey field 
dramatically affect the estimation of environment, with the impact of edge effects depending upon 
the adopted environment measure. These edge effects considerably limit the usefulness of smaller 
survey fields (e.g. the GOODS fields) for studies of galaxy environment. In even the poorest groups 
and clusters, redshift-space distortions limit the effectiveness of each environment statistic; measuring 
density in projection (e.g. using counts in a cylindrical aperture or a projected n th -nearest-neighbor 
distance measure) significantly improves the accuracy of measures in such over-dense environments. 
For the DEEP2 Galaxy Redshift Survey, we conclude that among the environment estimators tested 
the projected n th -nearest-neighbor distance measure provides the most accurate estimate of local 
galaxy density over a continuous and broad range of scales. 

Subject headings: methods:data analysis, methods:statistical, galaxies :high-redshift, galaxies:statistics, 
surveys, large-scale structure of universe 



1. INTRODUCTION 

The observed properties of galaxies have long been 
known to depend upon the environment in which they 
are located. For instance, red, non-starforming galax- 
ies (e.g. local ellipticals and lenticulars) are system- 
atically over-represented i n highly over-dense environ- 
ment s such as clusters (e. g.lDayis fc Geller|1976t|D resslerl 
ITflSnt IPo^tTna.u Oellerll1984t IBalogh et al.ll199SD . Re- 
cent studies have shown that the observed correlations 
between galaxy properties and environment are not lim- 
ited to the cores of rich clusters, but extend to less 
dense domains including the outer regions of clusters , 
galaxy groups, and the field (e.g. Balo gh et alJ Il999t 
Carlberg et alJ 120011 iBlanton et all l2003t iBalogh et all 



2004t ICroton et al.1 120051 . 



There are a variety of physical processes that can read- 
ily explain these observational trends, including the ac- 
tion of dynamical friction, tidal stripping, or gas pres- 
sure in dense environments. These mechanisms, in 
combina tion with the hierarchic a l model of galaxy for- 
mation llKauffmann et all 119931 iSomerville fc Primackl 
1999; ICole et a,U 12000 L in which galaxies form in less 
dense environments and are then accreted into larger 
groups and clusters, are generally consistent with the cur- 
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rent observations. From the empirical evidence, however, 
it remains uncertain in what environment (s), by what 
mechanisms, and on what time-scales galaxies evolve 
from a field-like population to a cluster-like population. 
Still, the strong correlation of local galaxy density with 
galaxy properties over a broad range of environment does 
indicate that it plays an important role in galaxy forma- 
tion and evolution. 

To study galaxy properties spanning a broad and con- 
tinuous range of local environment requires a thorough 
census of the 3-dimensional galaxy distribution over a 
large volume. Such data sets are only collected via large, 
systematic redshift surveys. For the nearby galaxy popu- 
lation, wide-fie ld spectroscopic an d photomet ric surveys 
(e.g. 2dFGRS. IColless et al.ll2001l and SPSS. lYork . etall 
2000) have provided excellent data sets for studying 
galaxy environments ranging from voids to rich clus- 
ters. Recent results from these large surveys have found 
that galaxy environments correlate strongly with the col- 
ors, luminosities, and morphologies of local galaxies (e.g. 



Balo eh et alJ2 004: Ho gg et all2 003: Bla nton et al. 200J 



Hog g et al.H 2004). 



With the advent of new high-reds hift surveys (e.g. the 
VLT- VIMOS Deep Survey rVVPSV lTTFevre et ahfeoM 
20041 and the DEEP2 Galaxy Redshi ft Survey (DEEP2), 



Davis et al.l 120031 iFaber et alJl2005j) . studies of galaxy 



environment will be able to extend beyond the local uni- 
verse. Such deep, high-redshift surveys will provide a 
representative snap-shot of the galaxy population and 
corresponding local densities when the universe was half 
its present age, thereby permitting an investigation into 
the influence of environment upon galaxy evolution and 
formation. That is, extending environment studies to 
higher redshift will enable a determination of whether 
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the correlations among galaxy properties observed in the 
local universe are the result of physical processes acting 
over the entire lifetime of the galaxy or whether the cor- 
relations were established during the early formation of 
the galaxy. 

While redshift surveys have grown in scale and studies 
of galaxy environment have increased in prevalence, few 
published tests have detailed the degree to which envi- 
ronment measures are affected by survey limitations. For 
instance, the confined sky coverage of surveys introduces 
geometric distortions (or edge effects) which bias local 
density measures near boundaries or holes in the survey 
field. Environment statistics can also be impacted by the 
redshift precision and target selection requirements of a 
given survey. Mock galaxy catalogs provide an excellent 
means for testing the biases introduced to a given density 
measure by effects such as these. 

In this paper, we test the applicability of several pop- 
ular density estimators within dee p redshift surveys uti- 
lizing the mock galaxy catalogs of lYan et alJ (|2004f) . In 
particular, we investigate the effects of redshift preci- 
sion, survey field edges, redshift-space distortions, and 
target selection. We also devote specific attention to the 
DEEP2 survey with the goal of identifying the optimal 
density measure for use within the survey. Within this 
study, we do not consider global measures of environment 
trends (such as correlation functions), but instead focus 
on measures which can estimate environmental proper- 
ties of individual objects. The outline of the paper is as 
follows. In the next section, we present the mock galaxy 
catalogs used to test environment measures at high red- 
shift. Subsequently (§3), we describe the environment 
measures to be tested. In §4, we examine the signifi- 
cance of redshift precision in determining local galaxy 
densities. In §5, we conduct a detailed analysis of edge 
effects with respect to each environment parameter. In 
§6 and §7, we then address the influence of redshift-space 
distortions and target selection on the various environ- 
ment estimators. In §8, the roles of completeness and 
the survey selection function are discussed. Finally in §9 
and §10, we conclude with a summary of the applicabil- 
ity of each environment measure at high redshift and a 
discussion of the suitability of current deep surveys to 
measuring local galaxy densities. 

2. SIMULATING DEEP REDSHIFT SURVEYS 

Begin ning with the Cen ter for Astrophysics Redshift 
Survey l)Davis et alJ 119821) . large redshift surveys have 
played a major role in studying galaxy properties, mea- 
suring cosmological parameters, and studying the large- 
scale structure of the universe. With improvements to as- 
tronomical instruments, local redshift surveys have bal- 
looned in size and surveys at high redshift [z ~ 1) have 
become possible. At present, high-z surveys take two 
forms: (a) obtaining precise (a z < 0.001) redshifts using 
spectroscopic observations of galaxies (e.g. DEEP2 and 
VVDS) and (b) using deep photometry in many pass- 
bands to make less precise (a z > 0.05) photometr ic red- 
shift estimates (e.g. COMBO-1 f. lWolf et al!200l . Each 
of these has its advantages and disadvantages. 

Even utilizing highly- multiple xed, multi-o bject spec- 
trographs (e.g. DEIMOS, iFaber et alJ 120031) on large- 
aperture telescopes, a deep spectroscopic redshift survey 
requires a vast amount of telescope time and is invariably 



limited in the number of galaxies for which it can measure 
redshifts. Slit or fiber collisions constrain the number of 
objects able to be targeted during a given exposure while 
the forest of OH sky lines in the optical and infrared plus 
instrument defects and limited signal-to-noise cause red- 
shifts to be missed for some percentage of targeted ob- 
jects. Spectroscopic surveys benefit from a higher level 
of redshift precision which permits studies of kinematics 
within galaxies and galaxy groups, while also measur- 
ing spectral properties such as emission-line equivalent 
widths. On the other hand, using an imager with a large 
field-of-view and observing in many passbands, less pre- 
cise photometric redshifts can be obtained for nearly all 
galaxies above a given magnitude limit in the targeted 
field. For this reason, photometric surveys are often able 
to build larger samples and are ideal for measuring the 
galaxy luminosity function or galaxy-galaxy lensing for 
which a high level of velocity accuracy is not necessary. 

In t his paper, we em ploy the simulated galaxy cata- 
logs of lYan et alJ ((2004) to model both photometric and 
spectroscopic surveys at z ~ 1. The simulations and 
all work in this paper employ a ACDM cosmology with 
= 0.3, Qa = 0.7, h = 1, and er 8 = 0.9. The mock cat- 
alogs are derived from N-body simulations by populating 
dark matter halos with galaxies accor ding to a halo oc- 
cupation distribut ion (HOD) function ijPeacock fc Smith 
120001 tScliak 2000) which describes the probability distri- 
bution of the number of galaxies in a halo as a function of 
the host halo mass. The luminosities of galaxies are then 
assigned according to the conditiona l lumin osity func- 
tion (CLF) introduced bv lYang et alJ l)2003f) . which de- 
scribes the luminosity function in halos of mass M. Mod- 
els for the HOD and the CL F are constrained from the 
2dFGRS luminosity function llMadgwick et al.l2002D an d 
two-point correlation function ( Madgwi ck et al J 12003) . 
By assuming that the manner in which dark matter halos 
are popul ated with galaxie s does not evolve from z ~ 1 
to z ~ (jYan et alJl2003|> . the mock catalogs are built 
using simulation outputs at corresponding redshifts. The 
simulated galaxy catalogs show excellent agreement with 
the lower-z (0.7 < z < 0.9) DEEP2 correlation function 
(iCoil et al.ll2004l) and the COMBO- 17 luminosity func- 
tion (Wolf ct all l2003fk they therefore should provide a 
realistic data set for studying measures of the environ- 
ment of galaxies at z ~ 1. For further details regarding 
the co nstruction of the mock catalogs, refer to lYan et alJ 

ESS). 

From the volume-limited mock catalogs, we are able 
to mimic a typical photometric redshift survey by se- 
lecting all galaxies above a given magnitude limit and 
applying to each galaxy redshift a random offset drawn 
from a Gaussian distribution with standard deviation a z . 
We utilize the DEEP2 survey as a model high-redshift 
spectroscopic survey. The volume-limited mock cata- 
logs are selected according to the DEEP2 magnitude 
limit of Rab < 24.1 and passed through the DEEP2 
target-selection and slitmask-making code, which is able 
to place approximately 60% of available t argets on slit- 
masks for spectroscopy ijDavis et al.ll2004|) . Finally, 30% 
of objects are randomly rejected to reflect a conservative 
redshift success rate of ~ 70%. The 12 mock catalogs 
cover fields of 120' x 30' in area with a total of ~ 120 
DEEP2 slitmasks tiling the 1 square degree. To simulate 
larger survey fields, we tiled multiple mock catalogs with- 
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TABLE 1 

SUBSAMPLES SELECTED FROM 
THE MOCK GALAXY CATALOGS 



Sample 


Description 


volume-limited 


full mock galaxy catalog; L > 0.1L* 


magnitude-limited 


R AB < 24.1 


DEEP2-selected 


R AB < 24.1; 




passed DEEP2 target-selection and 




slitmask-making criteria; random ~ 70% 




redshift success rate; high redshift 




precision (<r z = 0.0001) 



Note. — We present a list of commo nly used sub-samp les 
drawn from the mock galaxy catalogs of Van et al, 1 2004 ) . 



out overlap or discontinuity. Such large-field mocks were 
essential for studying edge-effects (§5) and for building 
large sample sizes. 

In each mock catalog, we have a complete tally of the 
total galaxy distribution down to a luminosity of 0.1L*, 
along with subsets of objects which pass the DEEP2 
target-selection criteria, were placed on a slitmask for ob- 
servation, and yielded a successful redshift. Such a cen- 
sus enables detailed study of the survey selection function 
and the manner in which slitmask-making and target- 
selection affect environment statistics. Throughout this 
paper, we utilize several subsets drawn from the mock 
catalogs as described in Table ^ Note that for each 
mock galaxy, the simulations provide accurate positions 
in both real-space and redshift-space. 

3. ENVIRONMENT MEASURES 

The environment of a galaxy is typically defined in 
terms of the density of galaxies located in its imme- 
diate vicinity. However, a variety of density measures 
are often employed in estimating environment. For 
example, many previous analyses have focused on the 
identification of predefined groups or clusters of galax- 
ies, which can be contrasted to those galaxies not in- 
habiting these over-dense regions - that is, the field 
popu l ation ( e.g. |Kuntschner et alJ l2002t Ivan den Berghl 
V2(M ILewis et 3,1.1120021 IChristieTrl EoOOy Another ap- 
proach is to instead derive a continuous measure of the 
galaxy density distribution, such as by measuring the 
distance to the n th -nearest n eighbor (e.g. iGomez et alJ 
12003 iMateus fc Sodrel 12004ft or by directly smoothing 
the observed gal a xy distribution on a fixed scale (e.g. 
Hogg et all 120031: iBeuing et aD 120021 iKauffmann et all 
2004J) . The underlying theme in each of these methods is 
that one requires a measure of the local number density 
of galaxies at the position of each galaxy in the sample. 

In our analysis, we focus on density estimators that 
do not rely on identifying galaxy groups or clusters in 
any way. Lumping galaxies into predefined classifications 
provides a poorly sampled range of galaxy environments 
especially when compared to continuous measures of en- 
vironment. At high redshift, where dense regions are 
commonly under-sampled and clusters and groups are 
less numerous, a more continuous definition of environ- 
ment is all the more desirable. Still, identifying galaxies 
in groups at hig h z is possible and has been tested in a 
separate paper l|Gerke et al.ll2005|) . In this analysis, we 



compare three popular density estimates: n} -nearest- 
neighbor distance, counts in an aperture, and the Voronoi 
volume. This set of measures is in no way presumed to 
be complete. Other promising methods for measuring 
local galaxy density, including using a Gaussian kernel 
to smooth the ga l axy distribution ove r a given scale (e.g. 
iHogg et alJl2003t iBalogh et al.ll200|) . are not discussed 
in this work. 

3.1. n th -Nearest-Neighbor Distance, D n and D p>n 

As first employed bv iDresslerl |l980), the local galaxy 
density can be estimated using the distance to the n th - 
nearest, spectroscopically-observed neighbor of a given 
galaxy. Often, redshift information is simply employed 
to exclude foreground and background sources - by re- 
stricting neighbors to a given velocity interval - and 
the nearest-neighbor distance is measured in projection. 
Commonly, the projected n th -nearest-neighbor distance, 
D Pt n, is expressed as a surface density, E n = n/ (irDp n ). 
Measuring nearest-neighbor distances in projection is 
particularly useful when stu dying the density of galaxies 
in gr oups and clusters (e.g. !Dressle:rlll980l [Lewis et all 
2002), where the appropriate velocity interval by which 
to exclude background and foreground galaxies may be 
selected according to the velocity dispersion of the group 
or cluster. In this manner, one may confidently exclude 
galaxies not associated with the group or cluster; further- 
more, as shown in §6, measuring in projection minimizes 
the impact of redshift-spacc distortions. 

For less dense environments or poorly sampled groups, 
there may be few neighbors within the selected veloc- 
ity interval, causing D p n to reflect the distance to other 
structures rather than the local density. In environments 
where working in projection is problematic, an alterna- 
tive is to compute the n th -nearest-neighbor distance in 
3-dimensions by searching in spherical apertures for the 
n th -nearest, spectroscopically-observed neighbor. Simi- 
lar to its projected counterpart, the 3-dimensional (3-d) 
n th -nearest-neighbor distance, D n , is often expressed as 
a number density, p n = (3n) / (An ■ Throughout this 
paper, all n th -nearest-neighbor distances are quoted in 
units of comoving h" 1 Mpc and the symbols D n and 
Dp jn are employed to denote the 3-dimensional and pro- 
jected n th -nearest-neighbor distances, respectively. 

To study the effectiveness of the 3-d and projected n th - 
nearest-neighbor distance measures at tracing the local 
density of galaxies in different environments, we compute 
both D n and D Pt7l for a DEEP2-selected sample consist- 
ing of 12,636 galaxies as drawn from a 120' x 60' mock 
catalog. In Figure ^ we compare the values of D n and 
D p>n for each galaxy in the sample as measured using the 
redshift-space galaxy positions to the measured value of 
D n as computed using the real-space positions for each 
galaxy, which should reflect the true local environment. 
We find that at high densities, where redshift-space dis- 
tortions are greater, the projected n th -nearest-neighbor 
distance is superior at tracing the real-space density of 
galaxies but still suffers greatly from peculiar velocities. 
On scales corresponding to intermediate- and low-density 
environments, the 3-d measure of D n is a slightly more 
accurate estimate of the true galaxy distribution. For 
a DEEP2-selected mock catalog, ~ 15% of the observed 
sample resides in the regime (log 10 (I?5) < 0.5) where 
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the 3-d n th -nearest-neighbor distance saturates and loses 
sensitivity. 
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typical sizes of groups in the mock catalogs (see Fig. 
|2J) . For these values of n, the sensitivity of the n th - 
nearest-neighbor measure is rather independent of n. In 
under-dense environments, the dependence on n is very 
weak, whil e in dense environments the strong clustering 
of groups iTPadilla et al.ll2004l ICoil et aljEfflofo arranges 
to curb D n for n > n gIoup . Also, in galaxy groups, 
rcdshift-space distortions are a much greater source of 
error in D n than small variations in the choice of n. 

For the projected n th -nearest-neighbor distance mea- 
sure, we test the sensitivity of D p ^ n using line-of-sight ve- 
locity intervals ranging from ±750 km/s to ±2000 km/s. 
As shown in Table [21 using a larger velocity interval 
by which to exclude foreground and background sources 
increases the accuracy of D p n in dense environments 
but also sacrifices sensitivity at low densities. We find 
that for a DEEP2-selected sample, a velocity interval of 
±1000 km/s — ±1500 km/s is best suited for a broad 
range of environments. Compared to photometric red- 
shift errors (in the best datasets, a z ~ 6000 km/s), 
the sizes of the tested line-of-sight velocity windows are 
small. However, larger velocity intervals sacrifice sensi- 
tivity on small scales and provide poorer measures of the 
local density about each galaxy; a window large enough 
not to be dominated by photometric redshift errors is also 
large compared to the typical length-scales of large-scale 
structure (e.g. the correlation length and typical void 
sizes), and thus provides a poor measure of environment. 



Fig. 1. — (Top) We plot log 10 (-D5) computed in 3-d from the 
redshift-space distribution of galaxies in a 120' X 60', DEEP2- 
selected mock catalog versus the 3-d log 10 (D5) values measured 
using the real-space galaxy positions. (Bottom) We plot pro- 
jected log 1() (-J 7j-D Pi 5) measured from the redshift-space distribu- 
tion of galaxies in the same DEEP2-selected mock catalog versus 
the 3-d logjQ (-D5) values measured from the real-space galaxy posi- 
tions. For each galaxy, the projected n th -nearest-neighbor distance 
was restricted to neighbors within ±1000 km/s along the line-of- 
sight and the distance has been scaled by a factor of \/3/2 for 
comparison to the 3-d measure of distance. The sample plotted 
consists of 12,636 galaxies and is restricted to the redshift range 
0.7 < z < 1.4 and to galaxies more than 4 h~ 1 comoving Mpc from 
the nearest survey edge. For each sample, the Spearman ranked 
correlation coefficient, p, is given for galaxies in the high density 
(log 10 (-D5) < 0.5) and the low density (log 10 (D$) > 0.5) regimes 
as defined by the real-space D5 values where the high density do- 
main is selected to include the most dense third of the sample. The 
dashed line in each plot follows a correlation of p = 1. 

In addition to a lack of sensitivity on given scales, 
the behavior of the projected and the 3-d n th -nearest- 
neighbor distances depends upon the choice of n. A 
measure of the n th -nearest-neighbor distance effectively 
smoothes the galaxy distribution in a non-linear fashion 
according to the adopted value of n. If n is chosen to 
be much larger than the richness of typical groups in 
the sample, then the n th -nearest-neighbor distances for 
galaxies in these groups will be pushed to erroneously 
high values, as it will reflect the distance to the next- 
nearest structure. In this work, we study both the 
projected and the 3-d methods for computing the n th - 
nearest-neighbor distance employing a variety of values 
for n. We limit most discussion to values of n = 2, 3, 5 
which correspond to th e sizes of small gro ups detected 
in the DEEP 2 survey ijGerke et alJ 12005(1 and to the 
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Group Richness 

Fig. 2. — We plot the distribution of group richness in a 
magnitude-limited (Rab < 24.1), 120' x30' mock catalog restricted 
to the redshift range 0.7 < z < 1.4 — open boxes. The crossed boxes 
give the distribution in group richness as observed in a DEEP2- 
selected sample. This observed group catalog is derived from the 
magnitude-limited mock according to the DEEP2 target-selection 
and slitmask-making processes described in §2. Each galaxy has 
been a ssigned to a group according to the prescription of Y an et all 
(2004). Group richness is given by the number of galaxies popu- 
lating the group or virialized dark matter halo. Since DEEP2 only 
samples the z ~ 1 galaxy population at a rate of ~ 45%, those 
groups with 7 or more members in the magnitude-limited sample 
are typically observed to have a group richness of 3 or more in the 
DEEP2 sample. 



3.2. Counts in an Aperture, C 

Another method for estimating the local galaxy den- 
sity is to cou nt galaxies within a fixed metric aperture. 
For example. iHogg et al.l l|2003|) count spectroscopically- 
observed galaxies in the SDSS within spheres of radii 
8 h^ 1 comoving Mpc centered on each spectroscopically- 
observed galaxy. In high-rcdshift surveys where the sur- 
vey field may cover ~ 1 square degree or less, such a large 
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TABLE 2 

Dependence of D v , n on velocity window 



n 


Ajj (km/s) 


Phigh 


Plow 


3 


750 


0.454 


0.786 


3 


1000 


0.540 


0.744 


3 


1500 


0.598 


0.659 


3 


2000 


0.600 


0.602 


5 


750 


0.384 


0.831 


5 


1000 


0.451 


0.808 


5 


1500 


0.496 


0.736 


5 


2000 


0.505 


0.673 



Note. — For a range of velocity intervals (Ad) from ±750 km/s 
to ±2000 km/s, we compute the projected n th -nearest-neighbor 
distance for the galaxies in a 120' X 60' DEEP2-selected mock 
catalog where Av is the velocity window over which to exclude 
foreground and background interlopers. Each set of D PiTl values 
are correlated against the real-space D n distances computed for 
the same mock catalog, as illustrated in Figure We present 
the Spearman ranked correlation coefficients, p, in the low density 
(Plow) an d high density (phigh) regimes. For n = 5, the division be- 
tween the low and high density regimes is made at log 10 (-D5) = 0.5, 
while for n = 3 the distinction is drawn at log 10 (-D3) = 0.4. 



spherical aperture will be dramatically affected by the 
survey edges (see §5). For instance, within a 30" x 30" 
field (i.e. 20 hr 1 comoving Mpc on a side at z ~ 1), 81% 
of spherical apertures with a radius of 1 hr x comoving 
Mpc will fit within the surveyed field, while for apertures 
of radius 3 h^ 1 and 5 h" 1 comoving Mpc only 49% and 
25% of the field, respectively, will be unaffected by edges. 
Furthermore, local studies indicate that larger apertures 
do not provide any advantage or additional information 
worth this high price. Both observational and theoretical 
studies, suggest that galaxy properties are more closely 
related to dark matter halo mass and small-scale envi- 
ronme nt than the large-scale environment of the gala xy 
(e.g. iLemson fc Kauffmann 1999; Blant on et al.l l2004). 

While choosing smaller spherical apertures would re- 
duce the amount of survey volume affected by edges, 
apertures smaller than approximately ±1000 km/s along 
the line-of-sight will not be sensitive to galaxy groups 
or clusters. The counts in an aperture measure effec- 
tively smooths the data on some adopted scale, thereby 
losing sensitivity on smaller and larger scales. In our 
analysis, we employ a series of cylindrical apertures mea- 
suring 1-2 h^ 1 comoving Mpc transverse (radius) and 
±500 km/s - ±2000 km/s along the line-of-sight. The di- 
mensions of our cylindrical apertures are chose n to match 
the t ypical sizes of halos in the simulations (|Yan et al.l 
2004). 

3.3. The Voronoi Volume, V 

The Voronoi volume is a geometric measure that 
has s een use from engine e ring and biology to astron- 
omy l|Ramella et all 120011 iMarinoni et alJl2002|K Un- 
like counts in an aperture, the Voronoi volume does 
not smooth the galaxy distribution in any way. It pro- 
vides a continuous, adaptive measure of galaxy density 
on all scales by measuring a unique volume about each 
spectroscopically-observed galaxy. 

As illustrated in Figure [21 the Voronoi partition of 
space is the three-dimensional analogue of the two- 
dimensional Dirichlet tessellation, in which a plane con- 



taining a set of data points is divided into a set of poly- 
gons, each containing one of the points. A Voronoi poly- 
hedron is the unique three-dimensional convex region of 
space surrounding a data point (the seed), such that 
within the polyhedron every point is closer to the seed 
than to any other data point. The faces of the Voronoi 
polyhedron are defined by the perpendicular bisecting 
planes of the vectors connecting the seed to its neighbors, 
where a seed's neighbors are those points connected to it 
by the Delaunay complex - the set of tetrahedra whose 
vertices are at the data points and whose unique, cir- 
cumscribing spheres contain no other data points. The 
Voronoi partition and Delaunay complex are thus geo- 
metrical duals of one another. 

Computing the Voronoi partition for a galaxy redshift 
survey provides a natural way to measure the local den- 
sity of galaxies, since the volume of a galaxy's Voronoi 
polyhedron will vary inversely with the distance to its 
closest neighbors. For this reason, the Voronoi volume 
associated with each galaxy serves as a natural param- 
eterization of that galaxy's environment. Galaxies in 
dense regions will have small Voronoi volumes, while iso- 
lated galaxies will have larger volumes. Our methods for 
co mputing the Voronoi partition are identical to those 
of lMarinoni et a I1 l|2002T) . and wc refer the reader to that 
work for computational details and for further discussion 
of the usefulness and historical context of the Voronoi 
partition and Delaunay complex. In this paper, we will 
employ the symbol V to denote the Voronoi volume of 
a given galaxy and all Voronoi volumes are measured in 
units of comoving (/i -1 Mpc) 3 . 




Fig. 3. — The two-dimensional Voronoi diagram (Dirichlet tessel- 
lation) for an array of points (that is, a distribution of galaxies) is 
illustrated. The points represent the galaxies and the open squares 
identify the Voronoi vertices. Note that the areas of the polygons 
depend strongly on the local density of points. Similarly, three- 
dimensional Voronoi cell volumes vary inversely with local density. 



4. REDSHIFT PRECISION AND TARGET SELECTION 
RATE: PHOTOMETRIC VERSUS SPECTROSCOPIC 
SURVEYS 

As discussed in §2, photometric and spectroscopic red- 
shift surveys differ in the precision with which they are 



G 



Cooper et al. 



able to measure galaxy redshifts. To test the significance 
of redshift precision in measuring local galaxy environ- 
ment, we have produced a variety of mock surveys with 
differing characteristics. First, we simulate two photo- 
metric redshift surveys which mimic the varying pre- 
cisions of the COMBO-17 photometric redshift survey 
quoted in the literature. Our first simulated photometric 
redshift survey adopts a magnitude limit of -Rab < 24.1 
and redshift uncertainty of a z « 0.02 which reaches 
equally deep and is more precise than the COMBO-17 
specifications of R VKna < 24, a z « 0.03 as given by 
iWolf et alJ (|2003^ . In addition, we simulate a photo- 
metric redshift survey with the same magnitude limit of 
-Rab < 24.1 and a redshift precision of a z » 0.05 as 
specified for COMBO-17 bv iTavlor etHI ^2003). Both 
magnitude-limited samples are drawn from the same 
volume-limited mock catalog and include 22,961 galax- 
ies covering a 120' x 30' field. Note that our assumed 
redshift uncertainties are lower limits to the redshift 
precision for the COMBO-17 survey. As discussed by 
iBell et aLl l)2004lh the photometric redshift precision for 
COMBO-17 depends strongly on the galaxy type and 
redshift; galaxies such as starbursts which lack a strong 
4000A-break yield redshifts with much greater uncertain- 
ties (a z ~ 0.1), while at higher redshifts K-correction 
uncertainties introduce systematic redshift errors. 

Running the same volume-limited galaxy catalog 
through the DEEP2 target-selection and slitmask- 
making software and assuming a conservative redshift 
success rate (see §2), we also produce a mock spectro- 
scopic sample with redshift precision of a z ~ 0.0001, 
mimicing the DEEP2 redshift survey. This DEEP2- 
selected spectroscopic sample includes 9,302 galaxies cov- 
ering the same 1 square degree field (sampling ~ 50% of 
galaxies to the magnitude limi t). To simulate the V VDS 
"deep survey" in the CDF-S (|Le Fevre et al]l2004j) . we 
randomly select 25% of objects to the same Rab < 24.1 
magnitude limit. This "VVDS-selected" sample is an op- 
timistic simul ation of VVDS. assu ming a 100% redshift 
success rate l)Vanzella et al.l 12004^1 and ignoring differ- 
ences in the bandpass used. Lastly, as a comparison sam- 
ple, we select the full magnitude-limited mock (22,961 
galaxies at Rab < 24.1) assigning redshifts according 
to the real-space positions of the galaxies as defined in 
the mock simulations. Each environment estimator (n th - 
nearest-neighbor distance, Voronoi volume, and counts 
in an aperture) is then computed on the photometric, 
spectroscopic, and real-space galaxy samples. 

For this comparison, we restrict our analysis to the 
redshift range 0.7 < z < 1.4 and to only those galaxies 
at transverse distances of greater than 4 h~ x comoving 
Mpc from the nearest edge in the survey volume. These 
restrictions make edge effects in both the redshift and 
transverse directions negligible, but do not introduce any 
selection biases (see §5). Note that the mock catalogs 
are not subject to interior edges; that is, the simulations 
cover a contiguous 1 square degree of sky with no holes. 

We find that the precision of even the best photometric 
redshifts is not sufficient to measure local galaxy environ- 
ments. Figure 0] shows the comparison between Voronoi 
volumes, V, as measured using the real-space galaxy po- 
sitions compared to those calculated using the observed 
redshifts for two representative surveys. Even assum- 
ing redshift errors as small as a z w 0.02, the environ- 



ment measured in a photometric redshift survey is insen- 
sitive for all but the very lowest density environments; 
the Spearman ranked correlation coefficient between the 
real-space and photometric measures of Voronoi volumes 
is p = 0.4. For the spectroscopic survey, redshift-space 
distortions introduce some scatter at high densities, but 
the overall distribution of environments is well measured. 
In all, the Voronoi volumes measured from the observed 
spectroscopic redshift distribution trace the real-space 
Voronoi volumes with much greater precision, yielding a 
correlation coefficient of p = 0.73. Very similar results 
are observed for the D n , D Ptn , and C environment esti- 
mators. 

In Tabled we expand our analysis to a better sampled 
range of redshift uncertainties. Even if the precision of 
photometric redshifts is greatly improved - by a factor of 
2 or 4 - we find that low-resolution spectroscopic surveys 
with galaxy sampling comparable to DEEP2 provide a 
significantly better trace of the 3-dimensional galaxy en- 
vironment. It is only at very high redshift precisions 
(a z < 0.005) and when measuring densities in projection 
that photometric redshift surveys are able to rival their 
spectroscopic counterparts as probes of galaxy environ- 
ment. 

Among the spectroscopic redshift surveys simulated, 
the greater sampling and improved redshift precision 
of the DEEP2 survey prove significantly superior to 
the VVDS in tracing the real-space density of galax- 
ies. On the other hand, at precisions better than the 
30 km/s uncertainty in DEEP2 redshifts, redshift-space 
distortions dominate the ability to measure local densi- 
ties and thereby limit any advantage of improved red- 
shift measurements (see Table 0J. Note that the galaxy 
samples in Table |3 transition from a magnitude- limited 
(-Rab < 24.1) sample at low redshift precision to mimic 
photometric or grism spectroscopic redshift surveys to 
a sample selected using the DEEP2 target-selection and 
slitmask-making procedures or a VVDS-like selection to 
simulate higher-resolution spectroscopic surveys, which 
have superior redshift precision but lower sampling rates. 

5. EDGE EFFECTS 

When measuring galaxy densities within any survey, 
one must always be careful of edge effects introduced by 
the limited area of sky covered in the survey. Even using 
the largest optical telescopes and an instrument with a 
generous field-of-view, a deep redshift survey is limited in 
its ability to cover large regions. Furthermore, to min- 
imize the effects of cosmic variance on the data set, a 
high-redshift survey is likely to spread the sky coverage 
over several fields. This limits the amount of contiguous 
sky coverage and increases the proportion of the survey 
area that is near an edge. In addition to the edges cre- 
ated by the chosen geometry of the survey field(s) , edges 
and holes can be created in the data set by effects such 
as bright stars in a field or problematic regions in photo- 
metric detection which prohibit any galaxies from being 
targeted there. To start, we will restrict our discussion to 
survey edges in the plane of the sky, but later discussion 
will address edges in the line-of-sight direction. 

The general effect of edges on each density estima- 
tor is to push the measurement towards lower density. 
To quantify the degree to which each environment mea- 
sure is affected by edges we compute each measure on 
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Fig. 4. — (Left) We plot a comparison of Voronoi volumes measured for the 4,539 galaxies in a 120' X 30', simulated DEEP2 spectroscopic 
redshift survey sample (cr z ~ 0.0001) which are more than 4 ft -1 comoving Mpc away from the nearest survey edge and within the redshift 
range 0.7 < z < 1.4. Plotted is the log of the Voronoi volumes as measured in the simulated DEEP2-selected sample versus the log of 
the Voronoi volumes as measured on the full magnitude-limited, real-space catalog. (Right) We plot a comparison of Voronoi volumes 
measured for the 9,310 galaxies in a high-precision (a z ~ 0.02) simulated photometric redshift survey sample which were more than 4 h~ 
comoving Mpc away from the nearest survey edge and within the redshift range 0.7 < z < 1.4. Plotted is the log of the Voronoi volumes 
as measured in the simulated photometric survey sample versus the log of the Voronoi volumes as measured on the full magnitude-limited, 
real-space catalog. In each plot, the Spearman ranked correlation coefficient, p, is given in the upper-left corner. The dashed lines follow 
a correlation of p = 1. While both samples arc influenced by redshift-spacc distortions (see §6) and the spectroscopic sample (left) suffers 
from poorer sampling, the inferior redshift precision of the photometric survey causes the line-of-sight galaxy distribution to smeared out 
irretrievably on all but the largest scales. 



TABLE 3 

Environment measurements in photometric and 
spectroscopic surveys 



Sample 


cr z 


PV 


PD 5 


P°P,5 


Rab < 24.1 


0.05 


0.307 


0.310 


0.389 


Rab < 24.1 


0.02 


0.400 


0.396 


0.473 


Rab < 24.1 


0.01 


0.494 


0.478 


0.575 


Rab < 24.1 


0.005 


0.596 


0.579 


0.688 


Rab < 24.1 


0.0025 


0.675 


0.677 


0.803 


VVDS-selected 


0.001 


0.640 


0.655 


0.689 


DEEP2-selected 


0.0025 


0.611 


0.625 


0.704 


DEEP2-selected 


0.001 


0.691 


0.716 


0.785 


DEEP2-selected 


0.0001 


0.726 


0.749 


0.802 


DEEP2-selected 





0.726 


0.751 


0.801 



Note. — For a range of redshift precisions (<r z ) we compute 
the 5 th -nearest-neighbor distance (-D5), Voronoi volume (V), and 
the projected 5 th -nearest-neighbor distance (D Pi s), for the galax- 
ies in a 120' X 30' mock catalog. As described in the main text, 
three samples are selected from the same mock catalog accord- 
ing to (a) the DEEP2 target-selection and slitmask-making pro- 
cedure [DEEP2-sclcctcd] , (b) a magnitude limit of Rab < 24.1, 
and (c) a 25% random sampling to the same magnitude limit of 
Rab 5: 24.1 [VVDS-sclcctcd] . Restricting to the redshift range 
0.7 < z < 1.4 and more than 4 hr~ 1 comoving Mpc removed from 
a field edge, the magnitude-limited sample contains 9,310 galax- 
ies while the DEEP2-sclcctcd and VVDS-sclcctcd samples include 
4,538 and 2,375 of those galaxies, respectively. For each sample 
and redshift precision, we (rank) correlate estimates of V, D5, and 
Dp, 5 with similar measures computed on a magnitude-limited sam- 
ple using the real-space galaxy positions. Note that u z = 0.0001 
corresponds to the redshift precision of the DEEP2 survey and that 
surveys reflecting the COMBO-17, VVDS, and DEEP2 attributes 
have been highlighted in bold font in the table. 



a large DEEP2-selected spectroscopic mock galaxy cat- 
alog covering a wide field, 120' x 90'. From the center 
of this larger simulation we extract a smaller rectangular 
survey field, covering 40' x 30', and rerun each environ- 
ment measure on this data set. For every galaxy in the 
smaller survey field, we then have measurements of en- 
vironment unaffected by edges (when measured on the 
larger sample) and measurements in which survey edges 
play a greater role (when measured on the smaller field). 
Trimming to the redshift range 0.7 < z < 1.4, the smaller 
sample consists of 2,803 galaxies. In the following sub- 
sections, we discuss how each environment measure is 
affected by the edges of the survey region on the plane of 
the sky. Note that we also ran tests incorporating holes 
and irregular survey edges with very little change in the 
relative results for the tested environment measures. 

For each galaxy in the studied sample, we express the 
difference in a given environment measure due to survey 
edges as a fraction of the width of the distribution for 
that measure. Specifically, we define the percent change 
in environment measure X for galaxy q by 

AeW= log 10 (^)-log 10 (X M ) 1Q0% 

01 

where ^2,9 is the measure of X for galaxy q computed on 
the smaller mock, X\ tq is similarly computed on the wide- 
field mock, and o\ is a measure of the Gaussian width 
of the logarithmic distribution of environment measure 
X as calculated in the larger simulation. Quantifying 
the change in each environment measure in this fashion 
enables the role of edge effects to be compared between 
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different environment estimators in a uniform manner. 

5.1. Survey Edges and n th -Nearest-Neighbor Distance 

The n th -nearest-neighbor distance environment mea- 
sure - in both projection and 3-dimensions - is affected 
by edges in a predictable manner. Any galaxy with 
an edge located closer than the measured n th -nearest- 
neighbor distance must be affected by the survey edges. 
However, to remove all such galaxies based on this cri- 
terion (D n > -Dodge or D PtU > £>cdgc) biases the sam- 
ple towards over-dense environments by excluding less- 
dense regions over a greater volume than more-dense re- 
gions. To avoid such biasing of the sample, a simple cut 
in edge distance can be made - excluding all galaxies 
within some distance of the nearest edge. This cut in- 
troduces no environment-dependent bias, but does allow 
some contamination to the sample at under-dense envi- 
ronments depending on the severity of the cut. In our 
simulations, we find that removing all galaxies within 
2 h^ 1 comoving Mpc of a survey edge creates a catalog 
with minimal contamination (roughly 5% of the sample 
has A e (D 3 ) > 10%) and still retains 65% of the data 
set. Relaxing the constraint to -D c dge > 1 h comoving 
Mpc, the level of contamination in the sample doubles to 
~ 10% with A e (£> 3 ) > 10% while the percentage of the 
sample retained increases to ~85%. 

As illustrated in Figure El we find that edge-effects 
show a clear dependence on n; the level of contamina- 
tion due to survey edges in the plane of the sky increases 
by roughly a factor of two for D§ relative to D3. Also, for 
a fixed value of n, the projected n th -nearest-neighbor dis- 
tance is slightly more robust to edge-effects in the regime 
where sample size is maximized (D < 2 h^ 1 Mpc in Fig- 
ure EJ). 

5.2. Survey Edges and Counts in an Aperture 

For an aperture of fixed comoving size, the edge effects 
upon the counts in an aperture density measurement are 
easily understood and cleaned from the sample. Only 
galaxies located within r t of an edge are affected, where 
2r t is the transverse diameter of the chosen aperture. 
Thus, by removing any galaxies within r t of an edge, the 
sample is entirely devoid of edge-affected galaxies. Such 
a trimming of the data set does not introduce a selection 
effect, that is, there is no bias towards environments of 
a given sort. 

In our simulated spectroscopic data set of 2,803 galax- 
ies, ~ 15% of the sample are positioned within r t of a 
survey edge using a cyclindrical aperture of r t = 1 hr x 
comoving Mpc. However, only for a scant ~ 3% of the 
sample did we find AC ^ when comparing measure- 
ments of C made on the smaller simulation to those made 
on the wide-field sample. One possible means for sal- 
vaging some edge-affected galaxies would be to scale the 
measured counts in each aperture by the amount of the 
aperture contained within the survey area. Due to the 
low rate at which C is perturbed by survey edges, how- 
ever, such a correction actually causes an overestimate of 
C for the majority of galaxies near the edge (D e d gc < ''t) 
of the survey field. 

5.3. Survey Edges and Voronoi Volumes 

Due to the geometrical complexity of the Voronoi tes- 
selation, understanding the effects of edges on the cal- 
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Fig. 5. — For the smaller sample of 2,803 galaxies, we plot the 
percentage of the simulated sample with -D c d ge > D {dotted line) 
where D e d gc is the projected distant from a given galaxy to the 
nearest survey edge. Note that the dotted line has been scaled 
down by a factor of 3 in order to fit within the plotting range. Of 
those galaxies meeting the edge-distance condition (-D c d gc > D), 
we also plot the percentage of galaxies corrupted by a survey edge 
on the plane of the sky according to the condition, A e (X) > 10%. 
For n = 3 and n = 5, the dashed lines and solid lines illustrate 
the contamination of D v , n and D n , respectively, as a function of 
the edge restriction Z5 e dge > D. In measuring D P: „, we employ a 
velocity interval of ±1000 km/s. 



culated Voronoi volumes is less straightforward than for 
the previously discussed density measures. For galaxies 
very close to exterior edges in the survey field, Voronoi 
volumes can be unbounded and such galaxies should be 
consequently discarded from the sample. On a more sub- 
tle level, edge effects - including interior edges - will also 
cause volumes to be increased in size while the volumes 
still remain bounded. Some of these edge-affected vol- 
umes can be detected as having Voronoi vertices outside 
of the survey field. However, many other edge-affected 
Voronoi volumes are not detectable in such a manner. 

As illustrated in Figure H3 even excluding galaxies lo- 
cated near a survey edge (e.g. within 2 h -1 comoving 
Mpc) , the distribution of Voronoi volumes is still greatly 
affected by edges with a bias towards large volumes be- 
ing pushed to even larger values (see Fig. EJ. It is pos- 
sible to minimize this effect by retaining only galax- 
ies with V below some limit. In our simulations, by 
truncating at log 10 (V^) = 3.1, the amount of contami- 
nation due to edge effects can be reduced to approxi- 
mately 20% of the sample with A e (V) > 10%. While 
making such cuts according to distance to the nearest 
edge (.Dodge > 2 h^ 1 comoving Mpc) and Voronoi vol- 
ume (\og 10 (V) < 3.1) effectively reduces the number of 
edge-affected galaxies in the sample, it also restricts the 
dynamic range of the Voronoi measure and considerably 
reduces the size of the sample. For our simulated spec- 
troscopic sample of 2,803 galaxies, the Voronoi volume 
density measure was the most dramatically affected by 
edges with >45% of the sample being corrupted, that is, 
having A e (V) > 10%. 
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Fig. 6. — (Top) For the 1,790 galaxies located more than 
2 h~ 1 comoving Mpc from an edge, we plot A e (V) versus the 
log of the Voronoi volume, V , as measured in the smaller field. 
Note that A e (V) is given by equation The clear trend is for 
galaxies in large volumes to be more severely affected by survey 
edges. The dashed line shows a proposed cut in Voronoi volume 
(log 10 (V) < 3.1) designed to remove highly edge-affected galaxies 
from the sample. (Bottom) We plot the distribution of Voronoi 
volumes for the full sample of 2,803 within the redshift window 
0.7 < z < 1.4 as measured on the wide-field mock (solid) and for 
the sample of 1,720 galaxies greater than 2 h~ 1 comoving Mpc of 
an edge and with log 10 (V) < 3.1 as measured on the smaller simu- 
lation (dotted). Note each histogram has been normalized to have 
an integral of unity. The overall distributions are in good agree- 
ment up to log 10 (V) = 2.7. The excess of galaxies with volumes 
2.7 < log 10 (V) < 3.1 as measured in the smaller simulation results 
from edge-affected galaxies contaminating the final sample. 



5.4. Effects of Finite Redshift Range 

As a secondary effect, the finite redshift range probed 
by any survey imposes edges in the linc-of-sight direction. 
The role of these edges is more easily handled by restrain- 
ing all scientific analysis to a limited, well-sampled range 
of redshifts. In the DEEP2 survey, the ability to measure 
rcdshifts at z > 1.4 or z < 0.7 decreases significantly as 
the [Oil] emission-line doublet leaves the observed op- 
tical window. In a DEEP2-selected mock catalog, we 
find that by restricting our sample to those galaxies at 
0.7 < z < 1.4, less than 1% of the sample has a 5 th - 
neareast-neighbor distance, D5, greater than the distance 
to the z — 0.7 or z = 1.4 edge. Similar contamination 
rates are found for the other environment estimators. 

A second concern for spectroscopic redshift surveys is 
the possibility of missing redshifts over specific wave- 
length intervals, especially in the far-optical and near- 
infrared where OH sky lines can dominate the spectrum. 
At z ~ 1 where optical surveys often rely upon a singular 
spectral feature (e.g. the [Oil] doublet at A rcs t ~ 3727A) 
for redshift measurements, a hole in wavelength cover- 



age translates directly into a hole in the survey's redshift 
sampling. For the DEEP2 survey, the high-resolution 
[R ~ 5000) of the DEIMOS data minimizes this effect, as 
the sky lines are then narrower than the components and 
spacing of the [Oil] doublet. In truth, the DEEP2 red- 
shift distribution exhibits no significant cross-correlation 
with a sky spectrum mapped to redshift according to 
either the central wavelength of the [ Oil] doublet or 
the w avelengths of either subcomponent ijNewman et alJ 
|2005[) . However, for lower resolution surveys such as the 
VVDS, windows of redshift insensitivity may be a con- 
cern that must be addressed in measuring galaxy densi- 
ties. 

6. REDSHIFT-SPACE DISTORTIONS 

While spectroscopic redshift surveys are able to mea- 
sure galaxy redshifts with great precision, redshift mea- 
surements by nature are measurements of velocity and 
not distance. Accordingly, converting differences in red- 
shift to relative line-of-sight distances is subject to the 
peculiar velocities of the galaxies. Such peculiar motions 
are greatest in dense regions such as groups or clusters 
where the velocity dispersion of the group causes the 
inter-member spacing to be larger in redshift space than 
in real space. Due to this environmental dependency of 
redshift-space distortions, it is essential to understand 
the manner in which they affect a given galaxy density 
measure. 

Within a mock DEEP2-selected spectroscopic galaxy 
catalog covering 120' x 90', we compute each environment 
measure using the both real-space positions of the galax- 
ies and the observed redshift-derived positions. Restrain- 
ing our analysis to galaxies at edge distances greater 
than Ah -1 comoving Mpc and within the redshift range 
0.7 < z < 1.4, we quantify the effect of redshift-space 
distortions on each environment estimator by calculat- 
ing the change in each measure as computed on the real- 
space mock relative to the corresponding measure de- 
rived from the observed spectroscopic mock. As in §5, 
we express the difference in a given environment measure 
as a fraction of the width of the real-space distribution 
for that environment measure. For example, the percent 
change in environment measure X for galaxy q is given 

by 

= log 10 -log 10 (X R , q ) 10Q% 

where X z q is the measure of X for galaxy q computed 
from the redshift-derived position and Ar i9 is similarly 
computed from the real-space position. The width, <jr, 
is determined via a Gaussian fit to the logarithmic distri- 
bution of environment measure X for all galaxies in the 
real-space simulation. Here, <7r can be measured on the 
real-space distribution of log 10 (X) or the redshift-space 
distribution with negligible difference for the DEEP2- 
selected sample. 

As illustrated in Figure the Voronoi volume and the 
3-dimensional 3 rd - or 5 th -nearest-neighbor distances are 
similarly affected by redshift-space distortions. For each 
measure, the effects of redshift-space distortions are non- 
negligible and as shown in Figure |H1 are greatest in over- 
dense environments. In comparison to the V and 3-d D n 
measures, the counts in an aperture density estimator, C, 
and projected n th -nearest-neighbor measure, D pn , are 
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less affected by the "fingers-of-god" due to their effective 
smoothing in the rcdshift direction; by definition, the 
projected estimators, C and D p n forfeit sensitivity in the 
redshift direction, which reduces their susceptability to 
rcdshift-space distortions. For nearly 80% of the sample, 
C is unaffected (AC = 0) by peculiar motions when using 
a cylindrical aperture with a length of ±1000 km/s and 
diameter of 1 h^ 1 comoving Mpc. The sensitivitity of 
C to redshift-space distortions is somewhat dependent 
upon the choice of the aperture size in the line-of-sight 
direction such that a smaller aperture is more adversely 
affected. For more than 90% of galaxies in our sample, 
we find | AC | < 1, again using an aperture with length of 
±1000 km/s. Similar results are found for the projected 
n th -nearest-neighbor distance measure; more than 80% 
of the sample meets the criterion A z (D Ptn ) < 5% for 
n = 3, 5 and using a velocity interval of ±1000 km/s to 
exclude foreground and background sources. 
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Fig. 7. — The effect of redshift-space distortions on the V and 
D n environment measures is illustrated. We plot the distribution 
of Aa(X) (see eq.gj for the 5,031 galaxies in the DEEP2-selected 
simulation within 0.7 < z < 1.4 and more than 4 h^ 1 comoving 
Mpc from a survey edge. The Voronoi volume and D n statistic 
are similarly affected by redshift-space distortions. Working in 
projection, for example using D Ptn or C, provides a much more 
robust estimate of galaxy density in over-dense environments. 



7. TARGET SELECTION AND OBSERVATION 

There are inevitable trade-offs between the number 
density of sources targeted for observation and the area 
of sky covered in a redshift survey. Clustering of high- 
redshift galaxies and fiber or slit collisions on multi- 
object spectrographs conspire to limit the fraction of 
target objects which a survey can observe at one time. 
Furthermore, not every object targeted will successfully 
yield a redshift, generally due to finite signal-to- noise and 
instrumental effects. 

The DEEP2 redshift survey will target ~50, 000 galax- 
ies covering 3.5 square degrees of sky, over 80 nights 
on th e Keck II telescope l|Davis et al 1 l2fMlFaEe7 et al.l 
2005). This impressive survey, however, will only target 
approximately 60% of available high-redshift (0.7 < z < 
1.4) galaxies in its four fields and successfully measure 
redshifts for about 75% of targeted galaxies. The DEEP2 
survey is designed with the goals of studying large-scale 
structure and galaxy properties at high redshift. Thus, 
the survey design attempts to maximize the number of 
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Fig. 8.— We plot A Z (D 5 ) (see eq. |2J as a function of Dr, as 
measured according to the real-space positions of galaxies. The 
sample plotted includes the 5,031 galaxies in the DEEP2-selected 
mock catalog with 0.7 < z < 1.4 and located more than 4 h~ 1 
comoving Mpc from a survey edge. Not surprisingly, redshift-space 
distortions are most problematic in over-dense environments. 



redshifts obtained, to sample the galaxy distribution over 
a broad range of length scales, and to minimize the effects 
of cosmic variance. Due to slit collisions on DEIMOS slit- 
masks, the DEEP2 survey systematically under-samplcs 
regions of sky with a high surface density of galaxies (see 
Figure^J) . It is critical for a study of galaxy environments 
to understand how this bias may affect the environment 
measured used. 

While the detailed effects of target selection and red- 
shift incompleteness are clearly specific to each survey, 
the goal of this section is to understand how a given envi- 
ronment measure is affected by the limited sampling com- 
mon to all deep redshift surveys. In this work, we adopt 
the DEEP2 survey as a representative high-redshift, 
spectroscopic survey. As discussed in §2, the DEEP2 
survey targets all galaxies at Rab < 24.1 according 
to a probabilistic algorithm which preferentially selects 
high-z galaxies. Applying the DEEP2 target-selection 
and slitmask-making algorithms to a magnitude-limited 
(-Rab < 24.1) mock catalog covering 40' x 30', the simu- 
lated survey targets and successfully measures redshifts 
for 2,839 of the 5,866 galaxies in the field and with 
0.7 < z < 1.4 (assuming a redshift success rate of ~ 70%). 
To study the combined effects of target selection, slit- 
mask making, and redshift success, we compute each en- 
vironment measure on the DEEP2-selected sample and 
on the full magnitude-limited sample. 

As illustrated in Figure ^| the counts in an aper- 
ture measure, C, shows no indication of an environment- 
dependent bias in the DEEP2 target selection. If DEEP2 
severely under-samples dense environments, then we 
would expect to see a saturation in the observed value of 
C at high densities relative to the estimation of C com- 
puted on the magnitude-limited sample. Instead, we find 
a linear relation extending to dense environments which 
follows the ~ 50% overall completeness of the DEEP2- 
selected mock catalog. 

Due to the fixed comoving aperture size of the counts 
in an aperture environment measure, it probes the same 
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Fig. 9. — {Top) We plot the distribution of angular-ncarcst- 
neighbor distances for galaxies in a 120' X 30', DEEP2-selected 
mock catalog. The solid line gives the distribution for all galaxies 
meeting the DEEP2 target-selection criteria. The dotted line plots 
the distribution for galaxies in the mock catalog which are selected 
for DEEP2 observation. (Bottom) The percentage of available tar- 
gets observed by DEEP2 as a function of angular-nearest-neighbor 
distance. The error bars are given according to Poisson statistics. 
The DEEP2 survey slightly under-samples dense regions of sky. 
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Fig. 10. — We plot the relationship between counts in an aper- 
ture, C, for the simulated DEEP2-selected redshift sample as mea- 
sured on the full magnitude-limited (R\b < 24.1) mock sample 
and as measured on the subset of the magnitude-limited sample, 
successfully observed by DEEP2. The contour levels plotted cor- 
respond to 1, 5, 25, and 100 galaxies. The dashed line follows the 
redshift rate of the DEEP2-selected mock sample (~ 50%). We 
employ a cylindrical aperture 2rt = 1 h _1 comoving Mpc in di- 
ameter and ±1000 km/s along the line-of-sight. All galaxies with 
D e( j ge < rt are excluded from the anaylsis so as to eliminate edge 
effects. No evidence is found for any environment-dependent bias 
in the DEEP2 target-selection and slitmask-making procedures. 



physical scale independent of how the targeted galaxies 
are selected. The n th -nearest-neighbor distance measure, 



on the other hand, can sample systematically different 
effective scales depending on the galaxy sampling. As il- 
lustrated in Figure im within the observed spectroscopic 
sample, the 5 th -nearest-neighbor distance is roughly trac- 
ing the 10 th -nearest-neighbor distance in the magnitude- 
limited mock; this is sensible, as the DEEP2-selected 
mock samples ~ 50% of galaxies and what was the 10 th - 
nearest neighbor in the magnitude-limited sample will 
typically be the 5 th in the DEEP2-selected sample. Sim- 
ilarly, the 2 nd - and 3 rd -nearest-neighbor distances are ef- 
fectively tracing the 4 th - and 6 th -nearest-neighbor dis- 
tances, respectively, in the magnitude-limited mock (see 
Figure ITT|) . While target selection and slitmask making 
affect the scale on which the n th -nearest-neighbor dis- 
tance samples the galaxy distribution, they do so in a 
manner which does not depend upon environment. Thus, 
while DEEP2 under-samples dense regions of sky, the 
survey does not under-sample dense environments (see 
Figure EH) . 

The limited galaxy sampling of the DEEP2 survey 
causes the calculated Voronoi volumes to be systemat- 
ically larger than if computed on the full magnitude- 
limited sample. We find that the level to which a given 
Voronoi volume is affected by the DEEP2 sample selec- 
tion does not depend on V or redshift; the limited sam- 
pling of the DEEP2 survey simply introduces a random 
scatter towards larger volumes. Similar to D n and C, we 
see no evidence for an environment-dependent bias due 
to the DEEP2 target-selection procedures. 

8. CORRECTING FOR THE SURVEY SELECTION 
FUNCTION 

For any magnitude-limited survey, the fraction of total 
galaxies in a volume-limited sample within the magni- 
tude limits varies - commonly decreasing - as a function 
of redshift. This variable sampling of the galaxy distri- 
bution as a function of redshift causes measurements of 
galaxy densities to depend strongly on z. For instance, 
if a survey under-samples at higher redshift, then esti- 
mates of D n and V at high z will be artificially inflated 
relative to estimates at low z; similarly C will be under- 
estimated at higher redshift. Often, magnitude-limited 
redshift surveys are trimmed to a volume-limited sub- 
sample to avoid these issues. At high redshift, however, 
this can dramatically reduce the sample size; for exam- 
ple, selecting a vol-limited sub-sample within a DEEP2- 
selected mock catalog excludes as much as 40% of the 
observed galaxies. Furthermore, over regimes where lu- 
minosity evolution is significant (A z > 0.1), even defin- 
ing a volume-limited sample can be problematic. 

To utilize the entire survey sample or for surveys that 
do not follow a simple magnitude- limited target selection, 
the variations in the galaxy sampling rate with redshift 
may be quantified in terms of a survey selection function, 
s(z), with which density estimates (number of galaxies 
per comoving volume or number of galaxies per projected 
comoving area) may be corrected as follows: 

X z (a, S, z) 



X (a,S,z) 



(3) 



s(z) ■ w(a, 6) ' 
where X z is the density estimate computed from the ob- 
served redshift distribution, w is a 2-dimensional survey 
completeness map which accounts for variation in red- 
shift completeness from field to field within the survey, 
and Xq is the corrected density estimate. 
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Fig. 11. — (Top) We plot the relation between D2 as mea- 
sured on the observed DEEP2-selected galaxy sample and D4 as 
measured on the corresponding magnitude-limited (Rab < 24.1) 
sample. (Middle) Plotted is the correlation between D3 as mea- 
sured on the observed DEEP2-selected galaxy sample and Dq as 
computed on the corresponding magnitude-limited (R\b < 24.1) 
sample. (Bottom) We plot D5 computed from the distribution of 
observed DEEP2-selected galaxies versus the D10 values measured 
on the magnitude-limited data set. In each plot, the Pearson corre- 
lation coefficient, r, is given in the upper left corner and the dashed 
line follows a correlation of r = 1. Again, no evidence is found for 
an environment-dependent bias in the DEEP2 target-selection and 
slitmask-making procedure. 



There are several ways to determine the selection func- 
tion ol a survey. The most common approach is to 
first estimate the galaxy luminosity function (LF) for 
all of the galaxies in the redshift survey, and then use 
it to predict the redsh ift distribution of the sample (e.g. 
Made wick et al.ll2003|h However, unless evolution is cor- 
rectly incorporated, the LF will not alone be able to cor- 
rectly predict the redshift dependence of the observed 
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Fig. 12. — (Top) We plot the distribution of "true" environment 
for DEEP2-selected galaxies in a 120' X 30' mock catalog where 
environment is traced by the 3-d 5 th -nearest-neighbor distance as 
computed within a volume-limited sample using real-space galaxy 
positons. The solid line gives the distribution for all galaxies meet- 
ing the DEEP2 target-selection criteria. The dotted line plots the 
distribution for galaxies in the mock catalog which are selected for 
DEEP2 observation. (Bottom) The percentage of available targets 
observed by DEEP2 as a function of environment. The error bars 
are given according to Poisson statistics. The DEEP2 survey does 
not under-sample dense environments. 



number density of galaxies in a deep survey. Further- 
more, working at high redshift, it becomes increasingly 
difficult to constrain the galaxy LF since observations are 
limited to the brightest sources, thereby making estima- 
tions of the characteristic luminosi ty, M*(z), and faint- 
end slope, a(z) : less secur e (e.g . iWillmer et all 120051: 
iWolf et all bona iBell et alJ 120041) . For this reason, the 
selection function for a survey is often estimated by 
smoothing the observed number density of galaxies, n(z), 
as a function of redshift and then normalizing according 
to an assumed dependenc e of the comovin g number den- 
sity of galaxies on z (e.g. lCoil et"al"ll2004j) . This has the 
disadvantage that density inhomogeneities in the survey 
will somewhat affect the derived redshift distribution, 
even with large smoothing kernels, due to the strength of 
cosmic variance; also, kernels large enough to minimize 
this will distort real features in n(z), especially where 
there are large gradients. On the other hand, any evo- 
lution in the observed number density of galaxies with z 
will be automatically incorporated into the estimation of 
s(z). 

Yet another approach to estimating s(z) is to compute 
an analytical fit to the observed dat a from which the se- 
lection function is then derived (e.g. iCoooer et al]l2005r 
iFaber et alJl2005|) . Similar to a selection function esti- 
mated from smoothing the observed n(z) distribution, 
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an analytical fit to the data - or "fitting" method for 
estimating s(z) - is subject to cosmic variance, but to a 
much smaller degree than the "smoothing" method, as 
small-scale variations in n{z) which do not match the 
model are not allowed. 

In this work, we estimate the survey selection func- 
tion for the mock DEEP2 survey according to four dif- 
ferent prescriptions: (a) estimating s{z) by smoothing 
the observed n(z) distribution in a DEEP2-selected mock 
catalog (120' x 30') assuming no evolution in the co- 
moving number densit y of galaxie s with redshift, using 
a similar algorithm as iCoil et all (|2004f> ( "smoothing" 
method), (b) fitting for the selection function assum- 
ing a functional form for the redshift dependence of the 
successfully-observed dN / dz and again assuming no evo- 
lution in the comoving number density of galaxies ("fit- 
ting" method), (c) determining the true selection func- 
tion by computing the number density of available tar- 
gets over many DEEP2 mock pointings relative to the 
volume-limited number density of galaxies in the mocks, 
and (d) deriving s(z) from the conditional LF assumed 
in constructing the mock catalogs. This last approach is 
identical to the commonly-used method of estimating the 
selection function using the measured LF and predicting 
the redshift distribution of the underlying galaxy popu- 
lation. The first two methods, (a) and (b), are analagous 
to methods one might use to derive s(z) solely from the 
observational data in a deep redshift survey and are sub- 
ject to cosmic variance and uncertainties in the assumed 
normalization and redshift dependence of the comoving 
number density of galaxies. The latter two methods are 
effectively identical, and test that the mock catalogs are 
working as advertised. 

In Figure ITTA we present the mock selection functions 
derived using each of the methods described above. In 
general, the agreement between the different approaches 
for determining s(z) is quite good. At the highest red- 
shifts (z > 1.0), the different estimations of the selection 
function differ due to differences in the assumed comov- 
ing number density of galaxies at high redshift. The 
conditional LF adopted in constructing the mock cat- 
alogs yields evolution that produces a decrease in the 
comoving number density, v(z), of galaxies at z > 1.0, 
while in estimating s(z) from the observed mock DEEP2- 
selected redshift distribution we assume a constant form 
for v{z). Both estimations of the comoving number den- 
sity at z > 1.0 are consiste nt with existing obse rvational 
evidence ijWolf et al.1120031: iWillmer et al]l2005h . 

The footprint of large-scale structure on a selection 
function derived from smoothing the observed DEEP2- 
selected mock n{z) distribution is reduced if a large 
smoothing kernel is used; here, we apply two successive 
smoothing windows of width Az = 0.15. If the smooth- 
ing kernel is too small, the presence or absence of struc- 
tures such as filaments or walls (i.e. cosmic variance) will 
cause us to overestimate or underestimate the fraction of 
galaxies sampled at a given redshift. Accordingly, over- 
or under-densities of galaxies will be inappropriately re- 
duced in amplitude when corrected by the survey selec- 
tion function; e.g. the presence of a filament will push the 
measured s(z) up at its redshift, reducing the corrected 
density measured, Xq, artificially. Any smoothing large 
enough to erase the effects of cosmic variance in a survey 
covering a few square degrees will, unfortunately, cause 
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Fig. 13. — {Top) We plot the galaxy redshift distributions for 
the mock catalogs. The solid line plots the observed redshift distri- 
bution, N(z), for galaxies within a DEEP2-selected mock catalog 
(120' X 30'). The dotted line gives the volume-limited distribution 
of redshifts for all galaxies (L > 0.1L») in the same 120' X 30' field. 
A probablistic roll-off in number density has been introduced into 
the mocks at z < 0.8 to simulate the high-z color-selection used in 
the DEEP2 target-selection. The dashed line shows the magnitude- 
limited (i?AB < 24.1) redshift distribution drawn from the same 
1 square degree field. (Bottom) The solid line gives the selection 
function, s(z), derived from successively smoothing the observed 
n(z) distribution in the DEEP2-selected mock DEEP2 catalog on 
scales of Az = 0.15 and assuming that ~45% of all available targets 
are successfully observed. The dashed line shows the analytical fit 
for the selection function as derived from the same DEEP2-selected 
mock. Lastly, the dotted line gives the estimation of s(z) derived 
from the conditional LF employed in constructing the mock cata- 
logs. Note that all of the s(z) curves have been normalized to peak 
at s(z) = 0.45. 



flattening in the shape of s(z), especially near the limits 
of the redshift range probed. Due to the drawbacks of 
smoothing, fitting for the selection function as detailed 
above is often a superior method for estimating s(z) from 
the observed data, but it does require assumptions about 
the form of dN/dz, which smoothing does not. 

To study the effectiveness of correcting the measured 
galaxy densities by the factor of l/s(z) (see eq. [3J , we 
have computed the projected 7 th -nearest-neighbor sur- 
face density, £7, within a volume- limited mock catalog 
covering 120' x 30' of sky. We then compare this to 
the projected 3 rd -nearest-neighbor surface density, S3, 
for those galaxies successfully observed in the DEEP 2- 
selected sample. We then correct these "observed" £3 
values using each of the s(z) shown in Figure and 
also attempt an empirical correction. This correction is 
given by dividing each observed density value by the me- 
dian E3 for galaxies at that redshift where the median 
is computed in a bin of Az = 0.04. Correcting the mea- 
sured densities in this manner converts the S3 values into 
measures of over-density relative to the me dian density 
and is simi lar to the methods e mployed bv IHogg et alJ 
( 2003) and lBlanton et all l|2003j) . 

Figure ITU illustrates the effectiveness of each selection 
function at reproducing the redshift dependence of the 
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galaxy density distribution as measured in the volume- 
limited sample. Within redshift bins of Az — 0.02, we 
compute the difference in median density between the 
corrected S3 values and the median density, £7, of ob- 
jects in the volume-limited sample. While each of the 
methods for estimating the survey selection function is 
an improvement over the uncorrected density distribu- 
tion, an empirical correction (as described in the previ- 
ous paragraph) which removes all z-dependence in the 
observed density distribution is at least as effective as 
correcting using a selection function (see Table @J|. 
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Fig. 14. — We plot the difference in the median density measured 
in a DEEP2-selected mock catalog — using the various s(z) estima- 
tions to correct the observed S3 values — relative to the median 
density, £7, computed in the corresponding volume-limited sam- 
ple, as a function of z. Ideally, these curves should be flat. The 
median densities are computed in redshift bins of width Az = 0.02. 
At redshifts z < 0.8, the £7 values measured in the volume-limited 
sample decrease due to the fall-off in v{z) forced in the mock cata- 
logs. For this reason, comparisons of the various s(z) are not valid 
at low redshift. The empirical correction (labeled as z- dependence) 
is at least as effective as correcting using a selection function. 



9. DISCUSSION 

Every environment measure that we have considered 
has its advantages and disadvantages. The counts in an 
aperture measure, C, lacks sensitivity in low-density en- 
vironments and while not lacking in dynamic range, it 
provides a non-continuous (or quantized) measure of en- 
vironment, a particular disadvantage if the typical value 
of C is small. It is best suited for working in dense envi- 
ronments where C is more robust to redshift-space dis- 
tortions than other measures and for analyses in which 
one wishes to classify a sample into coarse density bins 
or classes. The counts in an aperture statistic is unique 
among the environment estimators tested in that it mea- 
sures the galaxy density on a clearly defined, fixed length 
scale. In contrast, the projected and 3-dimensional n th - 
nearest-neighbor distance measures probe the radius en- 
closing some total number of galaxies and are not direct 
density measures. The C parameter also provides a great 
advantage via its robustness to survey edge effects. 

Similar to the counts in an aperture statistic, the 
projected n th -nearest neighbor distance measure is well 



TABLE 4 

Effectiveness of selection functions 



s(z) 


rms 


P 


observed 


0.180 


0.576 


smoothing 


0.094 


0.608 


fitting 


0.083 


0.603 


mock CLF 


0.096 


0.608 


z-dependence 


0.049 


0.606 



Note. — For each determination of the selection function, 
we compute the Spearman ranked correlation coefficient, p, be- 
tween the corrected 3 rd -neareast-neighbor surface density £3/5(2) 
measured in a DEEP2-selected mock catalog and the 7 th -nearest- 
neighbor surface density measured in a volume-limited sample. We 
als o pr esent the RMS fluctuations for each relation plotted in Fig- 
ure 1141 over the redshift range 2 > 0.8. Each of the selection func- 
tions improves the correlation (p) between the measured environ- 
ment within the DEEP2-selected sample and the "true" environ- 
ment measured in the volume-limited sample; the empirical correc- 
tion (z- dependence) has the smallest residual scatter with redshift 
(partially by construction). 
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Fig. 15. — We plot the correlation between projected 3 -nearest- 
neighbor distance, £>p,3, and the counts in a cylindrical aperture, 
C, measured in a DEEP2-selected mock catalog (40' X 30'). Here, 
we restrict the plot to those galaxies removed from any survey edge 
by (-Dedgc > 2 hT 1 comoving Mpc). The contour levels plotted 
correspond to N = 1, 10, 25, 50 and were computed using a sliding 
box as shown in the lower-left corner of the plot. The counts in 
an aperture utilize a cylindrical aperture of scale 2rt = 1 h^ 1 
comoving Mpc in diameter and ±1000 km/s along the line-of-sight. 
Similarly, a velocity interval of ±1000 km/s is used in restricting 
foreground and background galaxies in the computation of Dp, 3. 
In under-dense regions, the C values are found to saturate, thereby 
limiting the utility of the statistic. 



suited for measuring density in groups and clusters. 
However, unlike C, the projected D p „ parameter pro- 
vides a continuous estimate of galaxy density extend- 
ing to under-dense environments where it still provides 
a reasonably accurate measure. While slightly more ro- 
bust to edges than its 3-dimensional counterpart, D n , the 
projected n th -nearest-neighbor distance is more prone to 
survey edge contamination than the counts in an aper- 
ture statistic. Figure shows the correlation between 
D Pi 3 and C as computed in a 40' x 30' simualted DEEP2 
pointing. The saturation of C in less-dense regions is 
striking and proves to be a significant drawback for a 
density estimator which is otherwise extremely robust to 
survey edges and redshift-space distortions. 
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The 3-dimensional n -nearest neighbor distance and 
Voronoi volume statistics are the best suited for mea- 
suring under-dense environments. In groups and clus- 
ters, however, these density estimators are significantly 
affected by redshift-space distortions, causing each mea- 
sure to become saturated. As illustarted in Figure^O f ar 
removed from survey edges, the Voronoi volume and D 5 
measures agree very well over all environments observed 
in the DEEP2-selected mock catalogs. However, for the 
simulated DEEP2 survey data, the n th -nearest-neighbor 
distance is much more robust to edge effects and is less 
expensive to calculate. 
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Fig. 16. — We plot the correlation between Voronoi volume, 
V, and 3-d 5 th -nearest-neighbor distance, D5, measured on a 
DEEP2-selected mock field (120' X 30'). Here, we restrict the plot 
to those galaxies far removed from any survey edges, (D e< j ge > 
5 h —1 comoving Mpc). Note that the D5 values have been cubed 
to facilitate comparison to the Voronoi volumes. The Spearman 
ranked correlation coefficient, p = 0.85, quantifies the strong agree- 
ment between the different density estimators. 

For studies of environment at high redshift, includ- 
ing analysis in the DEEP2 survey, we conclude that 
among the environment measures tested the projected 
n th -nearest-neighbor distance provides the most accurate 
estimate of local galaxy density over the broadest range 
of scales. For work in dense environments, the D p ^ n offers 
great robustness to redshift-space distortions and main- 
tains a reasonably high level of accuracy in under-dense 
environments. While D p n can be affected by survey 
edges, contamination from geometric distortions is easily 
understood and effectively minimized without dramati- 
cally reducing the galaxy sample. 

10. CONCLUSIONS 

We have studied the applicability of several galaxy- 
density estimators within deep redshi ft surveys at z ~ 1 
utilizing the mock galaxy catalogs of lYan et alJ (|2004fl . 
We conclude as follows: 

1. Photometric redshifts derived from multi-band 
photometry (er 2 > 0.02) are not suitable for mea- 
suring galaxy densities. Current photometric red- 
shift surveys such as COMBO-17 do not have the 
redshift precision needed to study environment at 
high redshift. While more costly to obtain, spec- 
troscopic redshifts are requisite to accurately probe 
the local galaxy environment in a large survey. 

2. With the exception of the counts in an aperture 
estimator, C, survey field edges are a major source 



of contamination for each environment measure 
tested. To reduce these edge effects without bias- 
ing the sample, all galaxies within some comoving 
distance (~ 1 — 2 h^ 1 comoving Mpc for DEEP2) 
of a transverse survey edge should be rejected. At 
z ~ 1, excluding all galaxies within 1 h~ x comoving 
Mpc of an edge removes roughly 0.05 degrees along 
each dimension of the survey field. For smaller 
high-redshift surv eys, such as C FRS l)Lillv et alJ 
Il995h or CNOC2 llYee et al.H2000lh edge effects in- 
troduce contamination to a considerable portion of 
the survey data set, thereby limiting the statisti- 
cal power of the samples. Likewise, for a surve y 
of the GOODS-North field ijGiavalisco et alJl200l|) . 
edge effects would bias density measurements over 
~75% of the field. Testing each environment mea- 
sure on a simulated DEEP2-selected mock sample 
(40' x 30'), the Voronoi volume is most severely af- 
fected by edges, with more than 2 times as much 
contamination from edge effects than D n or D PtTl . 
The counts in an aperture measure displays the 
best behavior near edges of a survey field, with 
a nearly negligible portion of the sample contami- 
nated in our simulations. 

3. Redshift-space distortions are a significant and 
fundamental roadblock to measuring accurate 
galaxy densities in over-dense environments. The 
n th -nearest-neighbor distance measured in 3- 
dimensions and the Voronoi volume are most 
greatly affected, while estimators such as projected 
n th -nearest-neighbor distance and counts in an 
aperture - which smooth the galaxy distribution 
along the line-of-sight - are less affected by the 
"fingers-of-god" . Still, it should be noted that less 
than 15% of a simulated Rab < 24.1 galaxy sam- 
ple occupies environments at which the V and D n 
statistics saturate due to redshift-space distortions. 

4. The target selection algorithm employed by a sur- 
vey could lead to environment-dependent biases in 
the observed galaxy sample. The DEEP 2 survey, 
which slightly under-samples dense regions of sky, 
is equally sensitive at high and low densities. That 
is, we find that the DEEP2 survey equally samples 
all environments at z ~ 1 (see Figure H"2")l . Also, 
we find that while the DEEP 2 survey samples only 
~ 50% of galaxies at z ~ 1, this uniform incom- 
pleteness simply introduces a random scatter in the 
measured environments and does not introduce an 
environment-dependence bias. 

5. In examining the evolution of galaxy environments 
as a function of redshift, the estimation of the sur- 
vey selection function plays a critical role. Uncer- 
tainties in the comoving number density of galaxies 
at high z make comparisons over large redshift in- 
tervals (Az ~ 0.5) problematic. Apart from such 
ambiguities, simple empirical corrections for densi- 
ties as a function of redshift are highly effective. 

6. For the DEEP2 Galaxy Redshift Survey, the pro- 
jected n th -nearest-neighbor distance provides the 
most accuracte estimate of local galaxy density 
over a continuous and broad range of scales. The 
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D P n statistic is reasonably robust to redshift space 
distortions and still effective at tracing galaxy en- 
vironments in under-dense regions. 

7. Among current data sets at high redshift, we find 
the DEEP2 Galaxy Redshift Survey provides the 
best opportunity for measuring accurate galaxy 
environments over a broad and continuous range 
of scales. The high sampling rate and excellent 
redshift-precision of DEEP2 enable environments 
to be measured in even the most over-dense regions 
and yield improved accuracy over other deep sur- 
veys. Furthermore, DEEP2's high-precision red- 
shifts and large survey area (3.5 square degrees) 
minimize the effects of edges in both the transverse 
and line-of-sight directions. 
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